c=[1 2 3 4 5 6 7 8 9];
n=length(c);

A=circulant(c);
[V,D]=eig(A)
F=fft(eye(n))
D=fft(c);
ifft(diag(D)*F)